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Abstract. We examine the SLOCC classification of the non- 
normalized pure states of four qubits obtained by F. Verstraete et 
al. in [31j . The rigorous proofs of their basic results are provided 
and necessary corrections implemented. We use Invariant Theory 
to solve the problem of equivalence of non-normalized pure states 
under SLOCC transformations of determinant 1 and qubit permu- 
qq , tations. As a byproduct, we produce a new set of generators for 

the invariants of the Weyl group of type F4. We complete the de- 
, termination of the tensor ranks of four-qubit pure states initiated 

by J.-L. Brylinski [3J. As a result we obtain a simple algorithm for 
, computing these ranks. We obtain also a very simple classification 

---^ ' of states of rank < 3. 

1. Introduction 

We use the methods of Linear Algebra and Invariant Theory to study 
the problem of classification of pure quantum states of four qubits. 
Although we use the terminology common to Quantum Physics, we do 
not assume the reader is familiar with it, and we shall provide necessary 
definitions or references. We do not need a precise definition of qubits. 
It suffices to say that a qubit is a mathematical model for the quantum 
analog of an ordinary computer bit. A basic ingredient of this model 
is a 2-dimensional complex Hilbert space (see [2"fi] I2"7j). 

We shall work with four qubits. The Hilbert space of the k-th qubit 
will be denoted by Tik = C 2 with an orthonormal basis {eo,ei}. The 
Hilbert space for the quantum system consisting of four qubits is the 
tensor product 

H = Hi®H2®Hz® Hi. 
Occasionally we shall use Dirac's bra-ket notation in abbreviated form, 
e.g., 

\ijkl) = ti <E> &j ® ek <8> t\. 
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If Ak is an invertible linear operator on Tik, then we refer to A\ ® A 2 ® 
A% ® A 4 as an invertible SLOCC operation (reversible stochastic local 
quantum operations assisted by classical communication). 

A normalized pure state is a unit vector if) £ Ti up to a phase factor. 
However we shall work mostly with non- normalized pure states, i.e., 
mostly with nonzero vectors of Ti and we refer to them simply as pure 
states. 

The classification of pure states of three qubits is now well-known for 
both the group of SLOCC operations [JO] and the group of local unitary 
transformations [3J. The SLOCC classification of the pure states of four 
qubits was obtained by Verstraete et al. in [21]. However their list has 
an error which has not been noticed so far: the family L a b s is equivalent 
to a subfamily of L a b C2 . See Remark 13.51 for more detailed comments. 
The need to redo this classification on a more rigorous basis is also 
shared by some physicists [TU [22] . 

The study of the tensor ranks of four-qubit pure states was initiated 
in a recent paper of J.-L. Brylinski [3J, where he proposed that these 
ranks can be used as an algebraic measure of entanglement. We recall 
that the tensor rank of a pure state ip G Ti is defined as the least 
number r of product states whose sum is ip. By a product state we 
mean a (non-normalized) pure state of the form V\ ® V2 <S) v% ® V4. 
It is worthy of mention that the problem of calculating tensor ranks 
may be potentially very challenging. The case of tensor products of 
three vector spaces, of arbitrary dimensions, is relevant to the theory 
of algebraic complexity in Computer Science and we refer the reader 
to [1] for the exploration of this topic. 

Our objective in this paper is threefold. First of all we shall reprove 
Theorems 1 and 2 of [3TJ and at the same time improve and correct their 
formulations. Second, we give a simple method to test whether two 
semisimple states (see section [3J for the definition) are equivalent under 
SLOCC operations of determinant one and qubit permutations. The 
case of arbitrary pure states is also discussed. Third, we shall present an 
algorithm for computing tensor ranks for arbitrary pure states of four 
qubits. Along the way we shall give a different and simple classification 
of tensors of rank < 3, we shall determine the Zariski closure of the 
tensors of rank < 2 (a question left open in [3_] ), and we shall construct 
a nice set of generators for the algebra of polynomial invariants of the 
Weyl group of type F 4 (see Appendix B). 

There are several groups operating on Ti that are important for this 
paper. The most important ones are SLi oc = SL 2 x SL 2 x SL 2 x SL 2 and 
SL* oc = SLioc • Sym 4 , where SL 2 = SL 2 (C) is the special linear group 
in two dimensions and Sym 4 is the symmetric group on four symbols 
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which acts by permuting the four qubits. Similarly define GLi oc and 
GL* oc by replacing the group SL 2 with GL 2 (C). When we refer to M n 
or M mn as the spaces ofnxn resp. m x n matrices it is assumed they 
are over C, otherwise we will explicitly write M n (M), etc. 

Since we are interested in classifying the orbits of SL* oc , we shall 
need some basic facts about the polynomial invariants of this group. 
We denote by A resp. A* the algebra of complex analytic polynomial 
functions on 7i which are invariant under the action of SLi oc resp. 
SLj" oc . It is well known that A is a polynomial algebra in four variables. 
Its generators are algebraically independent homogeneous polynomials 
H, L, M and D of degree 2,4,4 and 6 respectively. See the paper 
[19] where these generators are explicitly constructed. The quadratic 
invariant (first constructed by Cayley, see [19] ) has the following simple 
expression 

i 

#(V0 = (- 1 ) i+i+ft ^,i,fc,oV'i-i,i-i,i-fe,i- 

i ,j,k=0 

The quartic invariants L and M are defined in section [3] and the sextic 
invariant D in section |H 

In section [2] we prove the basic result, Theorem 12. 8[ which gives the 
classification of indecomposable orthogonal representations of a very 
simple quiver Q (see Figure [1]). The matrix version of this result is 
stated in Theorem 12.101 One can view this theorem as a complex 
analog of the real SVD decomposition. The paper [31] contains such 
a theorem, for the case of square matrices only, but the uniqueness 
assertion for their proposed normal form is not valid. 

In section [3] we recall some basic facts about A. We also recall some 
important facts from the theory of infinitesimal complex semisimple 
symmetric spaces and introduce the notion of semisimple and nilpotent 
states ij} G Ti.. The main results are stated in Theorems 13.41 and 13.61 
They provide the classification of O4 x 04-orbits in M 4 and SLj* oc -orbits 
in Tt, respectively. The representatives of the 4 x 4 -orbits, organized 
in 17 families, are given in Table [TJ Nine of the families are selected 
to obtain a set of representatives of the SL* oc -orbits in H. The explicit 
expressions for these families are listed in Table [7| in appendix A. 

In section H] we show that A* is also a polynomial algebra in four 
variables and exhibit its generators H, T, E, II. These generators have 
degrees 2, 6, 8, and 12 and appear in a paper of Schlafli in 1852. We 
then show that these generators can be used to separate semisimple 
orbits. 
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The tensors if) G Tt of rank < 3 can be described by elementary 
means. This is accomplished in section The results obtained here 
are used in an essential way in the development of our rank algorithm. 

In section [6] we examine the nine families giving representatives of 
SLj* oc -orbits in 7i. For each of them we compute the tensor ranks of 
all states if) in the family. Finally, in section [7| we present a simple 
algorithm which computes the tensor rank of arbitrary pure states if) G 

n. 

In section [H] we summarize our results and make some comments 
on the problem of equivalence of two states under the group of local 
unitary operations. From the SLj* oc -classification of the pure states of 
four qubits it is easy to derive the GL* oc -classification. For a different 
approach to the GLi oc and GL* oc -classifications see the recent posting 
[18J on the arXiv. 

We thank A. Osterloh for his comments regarding the examples in 
section HI 

2. a complex analog of the real singular value 

Decomposition 

Our classification of SL* oc -classes of pure states of four qubits is 
based on a result of Linear Algebra which we discuss in this section. 

Consider M m>n (M) where we assume m < n. By A T we denote the 
transpose of a matrix A. If we consider the usual action of the real 
orthogonal groups O(m) and O(n) on M m>n (R), i.e, 

(2.1) (S 1 ,S 2 )-A = S 1 AS 2 1 

where Si G O(m), S 2 G O(n) and A G M m>n (R), then the orbits 
for this action are classified by diagonal matrices X G M m n (R) with 
non- increasing and nonnegative diagonal entries (the singular value de- 
composition theorem). 

There is a complex version of this theorem where M m>n (M.) is replaced 
by M m>n , the complex m x n matrices, and O(m) and O(n) replaced 
by the unitary groups U(m) and U(n). We are interested in another 
complex version of the above theorem where instead of U(m) and U(n) 
we use the complex orthogonal groups O m = O m (C) and O n = O n (C). 
We recall that O n is the subgroup of GL n (C) consisting of all complex 
matrices X such that X T X = I n . In the case m — n, such a theorem 
appears in the recent paper [31]. As the proof presented there is not 
completely clear, and the statement of the theorem we feel can be 
improved upon, we shall offer a different approach in this section. 

It is convenient to use the language of quivers. We just need one 
very simple quiver Q which has two vertices, say 1 and 2, and a single 
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Figure 1. The quiver Q 

O KD 

1 2 

We are interested only in orthogonal representations of this quiver. For 
simplicity we shall refer to them simply as representations. 

Definition 2.1. A representation of the above quiver Q is a 5-tuple 
(Vi, V 2 , 0i, 02, A) where V± and V 2 are finite dimensional complex vector 
spaces equipped with nondegenerate symmetric bilinear forms 0i and 
02 respectively, and A is a linear map V\ — > V 2 . We think of V\ and V 2 
as being attached to vertices 1 and 2 respectively while A is attached 
to the directed edge. 

Example 2.2. The most basic example of such a representation is 
given by V% — C n , V 2 — C m , where 0i and 02 are the usual dot products 
on C n and C m , and the linear transformation A : C™ — > C m is identified 
with an m x n complex matrix. Henceforth whenever we refer to C fc as 
a vector space we will assume that the bilinear form is the usual dot 
product. 

Let us recall some basic facts about representations of quivers in our 
setting. We start with the definition of a homomorphism of orthogonal 
representations. 

Definition 2.3. A homomorphism 

(2.2) S : (Vi, V 2 , 0!, 02, A) - (V[, V 2 \ 0;, 0' 2 , A') 

is a pair of linear maps Si : V\ — > V{ and S 2 '■ V 2 —>■ V 2 such that 

(2.3) / 1 (5 , i^i,S'iWi) = 0i(«i,u>i),V«i,wi G Vi 

(2.4) <j)' 2 (S 2 v 2} S 2 w 2 ) = 4> 2 (v 2 ,w 2 ),\/v 2} w 2 G V 2 
and S 2 o A = A' o S%. 

If S*! and S 2 are isomorphisms we say that S is an isomorphism. If 
there exists an isomorphism between two representations, then we say 
that the representations are isomorphic or equivalent. For instance two 
representations 

(2.5) A : C n -> C m ; £ : C" -> C m 
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are isomorphic if and only if there exist Si G O n and S 2 G O m such 
that S 2 A = BSi i.e., B = S^AS^ 1 . This means that two matrices 
A, B G M m n belong to the same orbit of O m x O n , acting on M m n in the 
usual way, if and only if the two representations above are isomorphic. 
Clearly, every representation of Q is isomorphic to one of the type given 
in example 12.21 

One defines the direct sum of representations in the obvious way. 
A representation (Vi, V 2 , 0i, (fi 2 , A) is nonzero if V\ or V 2 is a nonzero 
space. A nonzero representation is said to be indecomposable if it is not 
a direct sum of two nonzero representations. The Krull-Schmidt the- 
orem is valid, i.e., every representation decomposes as a direct sum of 
indecomposable ones and these indecomposable summands are unique 
up to permutation and isomorphism. 

Hence the classification of all representations of our quiver Q re- 
duces to the description of its indecomposable representations (up to 
isomorphism) . 

The three simplest nonzero representations are the following. 

Example 2.4. The representation A : — > C 1 is indecomposable. The 
matrix of the linear transformation A is the unique 1x0 matrix. 

Example 2.5. The representation A : C 1 —>■ is indecomposable. Its 
matrix is the unique 0x1 matrix. 

Example 2.6. Every representation A : C 1 — > C 1 is given by mul- 
tiplication by a fixed a G C. Its matrix is the lxl matrix [a]. If 
a = this representation is isomorphic to the direct sum of the repre- 
sentations from examples 12.41 and 12. 51 If a 0, the representation is 
indecomposable. 

The n x n symmetrized Jordan block J* (a) is the sum of the band 
matrix having a's on the diagonal and l's on the sub and super- diagonal 
and the matrix with i's on the opposite super-diagonal and — i's on the 
opposite sub-diagonal. For example 

a 1 i 0" 
1 a 1 + i — i 

1 +i a l-i . 

1 1 — i a 1 
-i la 

We point out that if n is odd, then the representation of Q given 
by J* (0) is decomposable. We have already seen this in Example 12.61 
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when n = 1. As a less trivial example consider 

1 + i 
1 + i l-i 
l-i 

By permuting rows and columns, we obtain 



' 1+i 








l-i 











l + l 


l-i 



Hence this representation of Q is the direct sum of the representations 
given by the two nonzero blocks, which are indecomposable. 

Remark 2.7. These blocks can be replaced with [l i] T and [l i] 
respectively. Indeed, we have [l + i 1 — z] = [l i\P where 

3 + i 3 — i 
3 — i —3 — i 

is orthogonal. We may simplify the symmetrized Jordan blocks as well. 
Note that for example J 3 X (a) is similar to 

'a 1 0" 

1 a i 
i a 

and since they are symmetric they are orthogonally similar. Hence, 
this matrix is in the same O3 x 03-orbit as J 3 x (a). 

The blocks J x (a) can be replaced by another kind of symmetrized 
Jordan blocks which consist of 3-diagonal symmetric matrices. For the 
description of these blocks see the recent paper [9]. 

We shall use a single matrix, as described in Example 12.21 to denote 
a representation of the quiver Q. With this in mind we can now state 
the following important classification theorem. 

Theorem 2.8. The representatives of the isomorphism classes of in- 
decomposable (orthogonal) representations of the quiver Q are given by 
the following matrices: 

(1) J* (a) for n > 1 where a 7^ if n is odd. The two values ±a 
give the same isomorphism class. 

(2) The (m + 1) x m matrix, m > 0, formed by using even index 
columns and odd index rows of J 2 x m +i (0) ■ 

(3) The transpose of the previous indecomposable. 



^(0) = 




cS 
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Proof. Let us first show that the representations A given in (1-3) are 
indeed indecomposable. Note that if the representations A and B of 
Q are isomorphic, then the matrices A T A and B T B are similar. Con- 
sequently, the number of indecomposable direct summands of A is at 
most equal to the number of Jordan blocks of A 7 A. 

In case (1), A = J* (a). If a ^ then A T A = A 2 has just one Jordan 
block and so the representation given by A must be indecomposable. 
If a = then n = 2m is even and A 2 is similar to J m (0) © J m (0). We 
leave to the reader to show that A must be indecomposable. 

The cases (2) and (3) can be handled together. Let A = J* (0) 
with n = 2m + 1 odd. Then A T A = A 2 is similar to J m (0) © J m +i(0) 
and so the representation A has at most two indecomposable direct 
summands. But we have seen in the discussion preceding this theorem 
that it indeed is a direct sum of two representations. Hence these 
summands must be indecomposable. 

We show next that every indecomposable representation of Q is iso- 
morphic to one of the representations (1-3). 

Let A : V — > W be an indecomposable representation where V = C n 
and W = C m . We have that A T A and AA T are linear operators on V 
and W respectively. Let us apply the Fitting decomposition 

(2.6) V = V ® Vi, W = W Q © Wx, 

where Vq and V± are v4 T A-invariant subspaces, A T A is nilpotent on Vq 
and invertible on Vi, and similar properties hold for Wo, W\ and AA T . 
Then it is easy to show that V ± V 1 , W ± Wi, and A(Vj) C Wi and 
A T (Wi) C Vi for z = 0, 1. 

This means that the representation A : V — > is the direct sum of 
the representations : V, — > where Aj is the restriction of A for 
z = 0, 1. As our representation is assumed to be indecomposable, we 
have V = W = or V x = W x = 0. 

Case 1: Vo = Wo = 0. Then m = n and A and A T are ismorphisms. 
By [S] A is a product of an orthogonal matrix and a symmetric 
one. A symmetric matrix is orthogonally similar to the direct sum of 
symmetrized Jordan blocks [UJ . Consequently, we can write A = PBQ 
where P,Q G O ra and B is the direct sum of symmetrized Jordan blocks. 
There is only one block, i.e, B = J£(cx) by the indecomposability 
assumption. As B is invertible, we have a ^ 0. It is also required to 
show that two symmetrized Jordan blocks J* (a) and J* (/?) are in the 
same isomorphism class if and only if a = ±j3. This can be seen as 
follows. If J£(a) and J^{P) gi ve isomorphic representations then 



J*(P)=PJ*(a)Q 
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so J* (a) 2 and J*(/?) 2 are similar and thus they have the same eigen- 
values hence (3 2 = a 2 . Conversely, note that J^(0) and — J*(0) are 
similar and symmetric so they are orthogonally similar. So there exists 
P G O n such that 



As —I n G O n we see that J* (a) and J* (—a) are in the same O n x O n - 
orbit, and so they give equivalent orthogonal representations of Q. 

Case 2: V\ = W± = 0. In this case the matrices A T A and AA T are 
nilpotent. This case occurs naturally in the theory of (infinitesimal) 
semisimple complex symmetric spaces. We omit the proof and refer 
the reader to [25118]. □ 

Remark 2.9. The classification of indecomposables in the above the- 
orem can be deduced from the general results on representations of 
symmetric quivers. We refer the reader to the recent paper of Derk- 
sen and Weyman J7j, where this new type of quiver is introduced and 
their representations (including the orthogonal and symplectic ones) 
are studied. In order to apply their results, our quiver has to be mod- 
ified by adding an additional directed edge from the second to the 
first vertex. The involution a, required by the definition of symmetric 
quivers, fixes the vertices and interchanges the two arrows. 

We can reformulate Theorem 12.81 in terms of matrices. 

Theorem 2.10. Let O m x O n act on M m , n in the usual way. Then the 
block- diagonal matrices 

A = Ay © A 2 © • • • © A k G M m>n , 

where each Ai is one of the matrices listed in Theorem \2.8\ are repre- 
sentatives of O m x O n -orbits. These representatives are unique up to 
permutation of the Ai 's and sign changes mentioned in that theorem. 

Remark 2.11. This theorem should be compared with [31, Theorem 1]. 
The authors consider only the square case m = n. Contrary to their 
claim, the canonical forms given there are not unique up to permutation 
of the diagonal blocks because some of their blocks have the shape 




4 x (o). 



PJ n x (a)P 



al n - J n x (0) 




R 
R 2 



i 
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and are made up of two of our rectangular blocks, one of size (p+1) xp 
and the other q x (q + 1). In the formulation given in Theorem 12.101 
there are more possible ways of combining such blocks, which leads to 
non-uniqueness. They also failed to mention that and J* (—a) 

belong to the same O n x O n -orbit. 

3. Classification of 4 x 4 -orbits in M 4 and SL^-orbits 

in n 

Our main objective here is to apply Theorem 12. 101 to the problem of 
SLj* oc -classification of pure states of TL. We start with some preliminary 
results, mostly well known. We have mentioned that a product state 
ip is of the form V\ (g) V2 ® V3 ® v±. We say a state ip is factorizable if, 
after a permutation of qubits, it can be written as the product of two 
tensors ip = ipi <8> ip2- As in [31j, for a tensor 

i 

i,j,k,l=0 

we define the 4x4 matrix ip by using the pairs ij as the row index and 
the pairs kl as the column index (we order these pairs as 00, 01, 10, 
11). By permuting cyclically the indices jkl we obtain two more such 
matrices ip' an d ip" ■ As in [19J, we denote their determinants by 

L = A 1234 , M = A 1342 , N = A 1423 , 

respectively. It is easy to verify that L + M + N = 0. 

Let Sk be the set of tensors with rank less than or equal to k. Surpris- 
ingly, it may happen that Sk is not Zariski closed. We shall denote its 
Zariski closure by Sk- We need the following result proved by Brylinski 
0. 

Proposition 3.1. The maximum rank of a tensor ip e 7i is 4. The 

affine variety S3 is irreducible and is defined by the equations L = M = 
0. Hence S3 has dimension 14. 

We prove the analogous result for S 2 , which was alluded to in [3]. 

Proposition 3.2. The affine variety S2 C 7i is irreducible of dimen- 
sion 10. Its ideal is exactly the ideal generated by the forty- eight 3x3 
minors of the matrices ip, ip 1 and ip" . 

Proof. Let I be the ideal generated by the forty-eight minors and W C 
TC its zero set. We used Singular [13] to verify that / is a prime ideal 
and that dim W = 10. On the other hand a simple computation (using 
Maple [23]) shows that the GL loc -orbit O of |0000) + | 111 1> also has 
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dimension 10. Since O C S 2 C W, we have C S 2 C W. As O is 
dense in W, the proposition is proved. □ 

Consider the action of SL 4 x SL 4 on M 4 given by: 

((P,Q),R)^PRQ T . 

The image of SL 4 x SL 4 under this representation is contained in 
SL(M 4 ), the special linear group of the space M 4 . This image is usually 
written as SL 4 <g> SL 4 , which means that we have two copies of SL 4 with 
their centers identified (glued together). 

The action of SLi oc on 7i gives rise to an action on M 4 via the map 
if) — > if). Explicitly, this action is given by 

((A u A 2 , Aa, A 4 ), $) -> (A 1 ® A 2 )iP(A 3 ® A 4 ) T , 

where A-y ® A 2 is the usual tensor product of matrices. 

The image of the first two factors SL2 of SLi oc under the action on 
M 4 is contained in the first factor of SL 4 <g) SL 4 . It is well known that 
this image is isomorphic to SL 2 <S> SL 2 = S0 4 but is different from S0 4 . 
We need to conjugate this image to obtain S0 4 . Clearly, the matrix 
which performs this conjugation is not unique. For that purpose we 
use the unitary matrix 



T 





1 








1 


1 





i 


i 










-1 


1 







i 








—i 



which we borrow from [31]. A slightly different such matrix Q is given 
in Makhlin's paper [2Tj. Now define R = R^ = Tif)T\ where the 
superscript ' indicates the hermitian transpose. It is assumed that the 
particular if> is obvious from the context when it is not written in the 
subscript. Finally we set 

R 

^~ -R T 



R — R 



-- At ® A 2 <g> A 3 ® Ai\4 
S0 4 are given by P x 



corresponds to = 
= T(A 1 ® A 2 )T^ and 



If A k e SL 2 then \<f>) = 
P X R^P 2 where P u P 2 G 

P 2 = T(A 3 ® A A ) T Tl Hence there is a 1-to-l correspondence between 
the SL loc -orbits in 7i and S0 4 x S0 4 -orbits in M 4 . 

For the following facts the reader can consult chapter 38 of [30J, 
and in particular Proposition 38.6.8. The SLi oc -orbit of if> is closed 
(in the Zariski topology) iff the S0 4 x S0 4 -orbit of is closed. It is 
well known that this is the case iff the matrix R^ is semisimple (i.e. 
diagonalizable). In this case we shall also say that if) is semisimple. 
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The Zariski closure of the SLi oc -orbit of ip contains the zero vector iff 
the same is true for the S0 4 x S0 4 -orbit of R^. Furthermore, this is 
the case iff the matrix is nilpotent. In that case we shall also say 
that ip is nilpotent. A nilpotent SLi oc -orbit, say O, is conical, i.e., if 
ip G O then also \ip G O for all nonzero scalars A G C. Hence O is also 
a GL loc -orbit. 

For any ip G T~l, the characteristic polynomial of R^ is given by 

t 8 + 2Ht 6 + (H 2 + 2L + AM)t A + 2(HL + 2D)t 2 + L 2 , 

where it is understood that H is short for H(ijf) etc. We will define the 
invariant D in the next section. If ip G S3 then L = M = and we 
obtain 

(3.1) t 2 (t 6 + 2Ht 4 + H 2 t 2 + AD). 

The discriminant of the cubic s 3 + 2Hs 2 + H 2 s + AD is equal to 

(3.2) 16D(H 3 -27D). 

The conjugation by the diagonal matrix J 4 © (— J 4 ) induces an invo- 
lutory automorphism 6 of 8 and its Lie algebra jj = so 8 (C) consisting 
of the skew-symmetric matrices in M 8 . Let 6 and p be the eigenspaces 
of 6 in q with eigenvalues +1 and —1, respectively. These eigenspaces 
consist of the matrices having the form 



V 0" 




"0 * 


★ 


resp. 


★ 



all blocks being of size 4. The space t is in fact a subalgebra of g, the 
Lie algebra of the subgroup K — 4 x 4 of 8 . The space p is a 
X-module with the action 

(\Pi 

This is an example of an (infinitesimal) semisimple complex symmetric 
space. The following theorem is a special case of general results about 
such spaces (301 Lemma 38.7.14]: 

Theorem 3.3. Let <p,ip G 7i be semisimple states. The invariants H , 
L, M , D take the same values at (f> and ip iff R^ and R^ belong to the 
same S0 4 x SO^-orbit, i.e., (ft and ip belong to the same SL\ oc -orbit. 

The special case m = n = A of Theorem 12.101 plays an important role 
in the sequel. We now state this special case in more detail. 

Theorem 3.4. The 17 families of matrices R, listed in Tabled classify 
the 4 x 4 orbits on M 4 up to permutation of diagonal blocks of the 





R 






-R T 


- 



p^p^ 1 

-P 2 R T P^ 
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same size and replacing the parameters a, b, c, d by +a, +b, +c, +d 
respectively. 

Table 1. i?-matrix representatives of 4 x 4 -orbits 





a 












a 






1. 




b 








2. 




b 








c 










c 


+ i 


1 
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1 b-i 
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-i 1 
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1 1 + i -i 
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a + i 1 
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1 a — i 
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a — i 
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11. 



a + i 1 






1 
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1 a — i 




12. 
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1 + i 


—i 




1 i 
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1-i 


1 








—i 





1 





" 1 


1 i 


—i 




1 i 
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1 + i 1-i 
-i 1 



1 


14. 


1 + i 1-i 
-i 1 












1 i 
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15. 



1 1 + i -i 

i 1 -i 1 



16. 





" 1 i 




17. 




1 i 













It should be noted that the representatives given in Table [T] may 
contain blocks different from those listed in Theorem 12.8} some of them 
have been simplified using Remark 12.71 Table [2] describes the Jordan 
structure of the R^ matrices. The lxl Jordan blocks are given by 
listing their eigenvalues. The symbol J2(±ic) indicates two 2x2 Jordan 
blocks with eigenvalues ic and — ic respectively, etc. The family 1 
consists of semisimple elements. On the other hand none of the other 
families contains a semisimple element. The nilpotent O4 x 04-orbits 
are easy to identify by using Table [2J Just set all parameters (if any) 
to in each of the 17 families. 

Table [3] gives a correspondence between the families of orbits found 
in [31] and those that we have identified. More precisely, for each of 
the 9 families given in [31] by explicit expressions we have determined 
the corresponding matrices and their Jordan structure as well as 
the corresponding O4 x 04-family in our notation (see Tabled]). The 
appearance of the imaginary units i in the expressions for eigenvalues 
of Rip is due to the fact that this matrix is skew-symmetric while the 
matrix P used in [31] is symmetric. 

Remark 3.5. Verstraete et al. [31] state that they found only 12 O4 x 
04-families, while we found 17. This is probably due to the fact that 
their Theorem 1 is not correct as stated. Their family L a t, 3 is equivalent 
to the subfamily of L abC2 obtained by setting c = a. We believe that 
there are two misprints in the formula for L a b 3 : the two + signs, in 
the last line of this formula, should be replaced by — signs. After this 
change, the family L a b 3 is equivalent to our family 6 and we have a 
perfect correspondence between their nine SLj* oc -families and ours. 

Some groups of families become one family once we examine how they 
behave under permutations of qubits. That is, an SLi oc -orbit from one 
family may be taken to an orbit in another family by permuting qubits. 
After this consideration there are nine different groups of families as 
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found in [H]. They are {1}, {2}, {3, 4, 5}, {6}, {7, 8, 9}, {10, 11}, {12, 13}, 
{14,15} and {16,17}. 

Table 2. Jordan structure of 

1. ±ia, ±z&, ±ic, ±id 2. ±za, ±ife, J 2 (±ic) 

3. ±ia, ±ifo, 0, J 3 (0) 6. ±ia, J 3 (±i&) 

4. Same as 3 

5. J2(±ia), J2(±ib) 



7. ±ia, 0, J 5 (0) 

8. Same as 7 

9. J 4 (±za) 

12. 0, J 7 (0) 

13. Same as 12 



10. J 2 (±ia), 0, J 3 (0) 

11. Same as 10 



14. J 3 (0), J 5 (0) 

15. Same as 14 



16. 0, 0, J 3 (0), J 3 (0) 

17. Same as 16 



Theorem 3.6. The orbits o/SL* oc on Ti. are classified by the nine fam- 
ilies 1,2,3,6,9,10,12,14 and 16 listed in appendix A (Table\7ty. Their 
R-matrices are given in Table\J\ States belonging to two different fam- 
ilies (from this list of nine) are not equivalent under SL* oc -operations. 
However, within the same family, different values of the parameters 
may give states belonging to the same SL^ QC -orbit. 

Proof. Denote by Ri the i?-matrix of the i— th family as given in Table 
[D Assume that k G {3,7,10,12,14,16}. One can easily compute 
the new i?-matrix, R' k , which results by applying the permutation 
(1,4) (2, 3) of the four qubits. Then it is easy to see that after mul- 
tiplying the first row and the first column of R' k by —1, we obtain 
exactly the transpose of R^ (if k = 3 or 7 this step is redundant). By 
inspection of Table [T] we see that R k = Rk+i- This means that the 
fc-th and (k + l)-st family of O4 x 0,4-orbits fuse into a single family 
of SL* oc -orbits. 

We leave to the reader to verify that the family 5 resp. 9 fuses with 
the family 3 resp. 7 into a single SL* oc -family. □ 

Let us point out that the O4 x 04-orbits may be disconnected and 
that different connected components may behave differently under qubit 
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permutations. We shall illustrate this in the case of the families 14 and 
15. These families are in fact single O4 x 04-orbits which we denote 
as (Di4 and i5 , respectively. Each of them has two connected compo- 
nents: 

014 = C[ A U 0[{, 015 = 7 015 U IJ 15 . 

These facts and some others that we will use are explained in [8]. To 
be precise, we assume that the Roman superscripts I and II are chosen 
so that the representative of 0m resp. 0is given in Table [H belongs 
to 0{ 4 resp. 7 0i5- The left resp. right multiplication of R by an 
orthogonal matrix with determinant —1 has the effect of switching 
the first two resp. last two qubits. The former leaves 0f 4 and 0f| 
invariant and switches 7 0is and /7 0is, the latter switches 0{ 4 and 0[{ 
and leaves 7 0is and 7/ 0i5 invariant. Switching qubits 2 and 3 is a 
new feature, not discussed in [SJ. We claim that in this case it has 
the following effect: The components 0{ 4 and 7/ 0is get interchanged 
while the components 0{ 4 and 7 0is remain invariant. To carry out this 
verification, one cannot rely on the Jordan structure of the matrices R 
as they are the same for both 0i 4 and 0i5. However these orbits have 
different afe-diagrams which makes it possible to verify the claim. For 
the definition of afe-diagrams see [T71 [251 IE]- Since the transpositions 
(1, 2), (2, 3) and (3, 4) generate Sym 4 , it follows that 0i4U0i5 is indeed 
a single SLj* oc -orbit. 

The redundancies mentioned in Theorem 13 . 61 will be addressed in the 
next section. 

Table 3. Correspondence of families of orbits 
Family in [21] Jordan blocks of R^ Our family 



Gabcd 


iiia, iib, iic, iid 


1 


Labc2 


Ha, iib, Ji(idc) 


2 


La2 62 


J 2 (iia), J 2 (±ib) 


5 


L a b 3 


iia, iib, J 2 (iia) 


? 


L a4 


J^iia) 


9 




J 2 (±ia), 0, J 3 (0) 


10 


L o 5e3 


•Mo), Mo) 


15 




0, J 7 (0) 


13 


L 3+ i,0 3+i 


0, 0, J 3 (0), J 3 (0) 
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4. Criterion for SLf r 



-EQUIVALENCE 



In this section we give a criterion for testing the equivalence of two 
states 0, ip e H under SLj* oc -operations. In the case when and ip 
are semisimple the criterion is very easy to use: one just has to verify 
whether the four invariants H,r,T,,U (the latter three to be defined 
below) take the same values on <p and ijj. 

As mentioned in the introduction, the algebra A of complex analytic 
polynomial functions / : 7i — > C which are SLi oc -invariant, i.e., satisfy 



is isomorphic to a polynomial algebra in four generators. Explicit gen- 
erators, as constructed in [19], are H,L,M and another polynomial D. 
The definition of D is somewhat involved. 
For j,k e {0, 1} let 



i,l=0 

where x ,Xi,y ,yi are independent commuting indeterminates. The de- 
terminant 



is a biquadratic form in the two sets of variables {xo,xi} and {yo,Ui}- 
There is a unique 3x3 matrix B such that this form can be written as 



Then D{iji) = det B G A and it is homogeneous of degree 6. By 
permuting cyclically the last three indices of ?/>, we obtain two more 
such invariants which we denote by E and F. 

One can easily verify that the four homogeneous polynomials 



are algebraically independent and invariant under the action of SL* oc . 
The degrees of these polynomials are 2,6,8 and 12, respectively. These 
polynomials appear in a work of Schlafli in 1852, who also noticed their 
invariance property under permutations of indices |29j . 

Let us now examine the mentioned redundancies of Theorem 13.61 
The most interesting case is that of family 1, the family of all semisim- 
ple orbits. The question we raise is the following: When are two states 



/(0"0) = /(VO> VgeSUnMeH, 



i 




^00 #01 
#10 #11 




H, T = D + E + F, E = L 2 + M 2 + N 2 , 
11= (L — M)(M - N)(N - L), 
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ip abcd and ip a ' b ' c 'd' in the same SLf oc -orbit? (By ip abcd we denote the 
state whose i?-matrix is the first matrix in Table [TJ) 

Let a be the subspace of Tt consisting of tensors if) with a diagonal 
matrix. If we identify TC with p using the map ip — > then a is a 
maximal abelian subspace of p consisting of semisimple elements. Such 
a subspace is known as a Cartan subspace of p. We mention that 
all Cartan subspaces of p are conjugate by SO4 x SO4, the identity 
component of the subgroup K = O4 x O4 of 0$. 

Let Na resp. Za be the subgroup of SLi oc which leaves a globally 
resp. pointwise invariant. Define similarly the subgroups iV£ and Z a 
of SL* oc . The quotient groups W a = N a /Z a and W£ = N^/Z* a act on 
a effectively. 

Let us use the diagonal entries of R^ as coordinates in a. It is easy 
to see that Wa can permute arbitrarily the coordinates a,b,c,d and also 
replace them with ±a, ±6, ±c, ±d provided the number of " — " signs 
is even. We conclude that Wa has order at least 192. On the other 
hand, a is a Cartan subalgebra of so$ and Wa is a subgroup of the Weyl 
group of the pair (sos, a). Since so$ has Cartan type -D4, this Weyl group 
has exactly order 192. We conclude that Wa coincides with this Weyl 
group. 

It is easy to check that Sym 4 C N^, i.e., all qubit permutations map 
a into itself. It is also easy to check that the permutation (1,2) (3, 4) 
acts trivially on a. It follows that the Klein four-group V < Sym 4 also 
acts trivially. The transposition (2, 3) sends the point (a, b, c, d) to the 
point 

— (a + b + c + d, a + b — c — d, a — b + c — d, a — 6 — c + d), 

i.e., it acts as the reflection in the hyperplane a = b + c + d. By using 
GAP [12], one can easily check that Wa and this reflection generate a 
group of order 1152. Clearly this is the Weyl group of type F 4 . 

Lemma 4.1. W^ is the Weyl group of type F4, of order 1152 = 2 7 ■ 3 2 . 

Proof. We have 

(4.1) [N* a : Z a ] = [N* a : N a }[N a : Z a ] = [N* a : Z* a }[Z a : Z a ]. 

We have seen above that [N a : Z* a ] = \W^\ > 1152. Since the Klein 
four-group V < Sym 4 acts trivially on a, we have [Z a : Z a ] > 4. On 
the other hand, since N a = N a n SL loc , we have [N a : N a ] < [SL^ oc : 
SL loc ] = 24. Recall also that [N a : Z a ] = 192. The equality 14. II implies 
now that all these inequalities are in fact equalities. In particular, 

TO = [N* a : Z* a ] = 1152. 
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Since Wa is a finite irreducible reflection group of rank 4, it must be 



Theorem 4.2. The restriction homomorphism p : A* — > B from the 
algebra A* to the algebra B of polynomial invariants on a is an 
isomorphism of graded algebras. The algebra A* is generated by the 
four homogenous algebraically independent polynomials H , V , E, and 
IT of degree 2, 6, 8 and 12, respectively. 

Proof. If / G kerp, i.e., / G A* and / vanishes on o, then / vanishes on 
all semisimple elements of TC. But the semisimple elements are dense 
in TC, and so / = 0. This shows that p is injective. 

In view of Lemma 14.11 we can apply to some well known facts 
about finite reflection groups, see for example [MJ Section 3.7]. The 
algebra B is isomorphic to a polynomial algebra in four variables and 
it is generated by four homogeneous polynomials of degree 2, 6, 8 and 
12. Moreover any set of four homogeneous generators of B must have 
these degrees. Now recall that the SL* oc -invariants H, T, E, II have 
exactly these degrees. Since they are algebraically independent, and p 
is injective, their restrictions to a are also algebraically independent. 
As their degrees are 2, 6, 8 and 12, they must generate B. Hence p 
is also surjective. We can now prove the following analog of Theorem 



Theorem 4.3. Two semisimples states <j>,ip ETC are SL* oc - equivalent 
iff the invariants i7 7 T ? E and IT take the same values at <p and ip. For 
arbitrary states <p,ip G TC, if at least one of the invariants H , V, E ; II 
takes different values on <p and ip, then the SL* oc -orbits of <p and ip are 
different. 

Proof. The second assertion is obvious. For the first assertion, we need 
only prove that its condition is sufficient. Assume that the condition is 
satisfied. Let a C TC be the Cartan subspace introduced above. Since 
every semisimple element ip E TC is SLi oc -equivalent to an element of 
a, we may assume that <fi,ip G a. By Theorem 14.21 and our hypothesis, 
all invariants of take the same values on <fi and ip. Since is a 
finite reflection group, we conclude that (ft and if) are PV^-equivalent. 
Since = N^/ Z* a , it follows that 4> and ip are ^-equivalent. In 
particular, they are SL* oc -equivalent. □ 

We can now use this theorem to show in a straightforward manner 
that certain tensors are not SL* oc -equivalent. In some situations one 
may use certain Bell inequalities as in [33] to show that two states 




the Weyl group of type F 4 . 



□ 



□ 
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Table 4. SLj* oc -invariants of some pure 4-qubit states 



State 


H 


r 


s 


n 


GHZ 


i 

2 











W 














10) 








1 

128 


i 

2048 


10') 








1 

128 


1 

2048 


lx> 








1 

128 


-1 
2048 



are not equivalent. However we feel that simply calculating the SLj* oc - 
invariants is a more straightforward approach and the above theorem 
will be enough for the majority of situations. 

Let us look at some examples. The generalized GHZ and W states 
in four qubits are 

^(loooo) + |1111)) 

and 

(4.2) i(|0001> + |0010) + |0100) + |1000)) 

respectively. Two important states in quantum teleportation [33], Eq. 
(22) and Eq. (2)] are the cluster state 

10) = i(|/? + o/?+o> + |/3 + o/ri) + 1/3-1/3-0) + |/ri7? + i)) 

where j/?*) = ^(|0) ± |1)) and the state \x) given by 

1 

2^/2 



\x) = -^=(|0000) - |0011) - |0101) + |0110) 



+ |1001) + |1010) + |1100) + |1111>) . 
There is also another cluster state mentioned in |2] Eq. (4)]: 

\<P') = i(|0000> + |0011) + |1100) - |1111». 

In Table H] we tabulate the values of H, T, E and IT on these five 
states. It is clear from this table that all five states belong to different 
SLj" oc -orbits except possibly for the pair \<f>) and | </>'). These two states 
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share the same invariants and, as both are semisimple, they indeed 
belong to the same SL* oc -orbit. 

We remark that the states e l?r//4 |0) and \x) belong to the same SL* oc - 
orbit because they are both semisimple and share the same invariants. 
As A. Osterloh has pointed to us, from the quantum physics point of 
view, the unit vectors |0) and e J7r//4 |0) represent the same pure state 
because they differ only by a phase factor. On the other hand, they 
belong to different SL* oc -orbits. The apparent discrepancy is explained 
by the fact that we are classifying the nonzero vectors in TC rather than 
the genuine pure states (see the Introduction). 

We can now sketch our procedure that one can use to decide whether 
two arbitrary pure states 0, ip G TC are SLj* oc -equivalent. It is under- 
stood that equivalence will mean SLj* 0( ,-equivalence for the rest of this 
section. Clearly if and ip are equivalent and is semisimple or nilpo- 
tent then ip must have the same property. 

Step 1. We compute the values of H, T, £ and II at and ip. If 
they do not agree then and ip are not equivalent. From now on we 
assume that they do agree. If and ip are semisimple, they must be 
equivalent. We shall now assume they are not semisimple. 

Step 2. Assume and ip are nilpotent and compute the Jordan 
structures of and R^. By inspecting Tabled with all eigenvalues 
set to 0, we see that apart from one case the Jordan structure of the 
i?-matrix determines uniquely the SL* oc -orbit. The exceptional case is 
when the Jordan blocks are of size 1, 1, 3, 3. Then there are two orbits. 
They can be distinguished by using a6-diagrams. One of these orbits 
is that of the generalized W state and the other is in family 16 (or 17). 
From now on we assume that and ip are not nilpotent. 

Step 3. The families 5 and 9 can be distinguished from the famlies 
3, 4 and 7, 8 respectively by the sizes of Jordan blocks of i?-matrices. 
By permuting qubits in both and i/j we can assume that both states 
belong to one of the families 2, 5, 6, 9, 10 or 11. After this reduction 
and ip are equivalent iff they have the same Jordan structure. 

5. Classification of tensors of rank at most three 

Here we provide some normal forms for tensors of rank 1,2 and 3 
under the action of SLj* oc and investigate some of their properties. A 
rank 1 tensor is just a product state and so it is in the same orbit as 
|0000). 

For the rank 2 case we have a few more situations to consider. Let 

■0 = 0-1 ® 0-2 ® a 3 ® a 4 
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Figure 2. Patterns for rank 2 tensors 
(a) * • • • (b) ] • • • (c) ( ( * ' 

■ •••• 1 1 • ■ 

be a rank 2 tensor. We may consider where linear dependencies occur 
amongst the sets Since we also consider the action of Sym 4 

there are really only 3 cases as in Figure [2j A line connecting two 
points in the z-th column means those corresponding two vectors are 
scalar multiples of each other. In case (a) let g% G GL2 be such that 
gi(ca) = Vit§ and g%{bi) = z^ei where vi is chosen so that gi G SL 2 , for 
i = 1,2,3,4. Hence the tensor reduces to a(|0000) + 1 1 1 1 1 ) ) , where 
a = v\v<iVzV/±- If is linearly dependent then we may set = 1 

since we are free to choose how the transformation acts on some other 
vector that creates a basis. Let us also describe how to handle (b). 
Since {ai,&i} is linearly dependent, we may assume 61 = a\. Hence 
the tensor is in the same orbit as a|0) ® ( 1 000) + 1 1 1 1 ) ) . Similarly (c) 
is in the same orbit as a|00) <8> ( 1 00) + |11)). Note that conversely any 
tensor in one of these forms is of rank 2. 

Proposition 5.1. A rank 2 tensor ip G 7i is SL* oc -equivalent to one 
of the following: 

(a) a(|0000) + |1111)), a^0, 
{b) |0)®(|000) + |111)), 
(c) |00)® (|00)+ 11)). 

In case (a) ip is semisimple and non-factorizable, while in cases (b) and 
(c) it is nilpotent and factorizable. 

Proof. The first assertion follows from the above discussion. We can 
assume that a = 1 in cases (b) and (c) since they are nilpotent orbits 
and so SL* oc and GLj* oc -orbits coincide. The second assertion is easy to 
verify. □ 

The case of a rank 3 tensor is not as easy to breakdown. The com- 
plications arise because now we have 3 vectors a^b^Ci being mapped 
under a SL^-transformation but can only control where 2 of them are 
mapped to in most cases. Let 

ijj = ai ® a 2 <g> a 3 <g> a 4 
+ b 1 ®b 2 ®h®b i 
+ ci ® c 2 <E> c 3 ® c 4 



be a rank 3 tensor. 
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FIGURE 3. Patterns for rank 3 tensors 




We can bring ifj into a reduced form by using SL 2 transformation on 
each qubit. The first three qubits are handled a bit differently than 
the last one as we shall see. We outline how to construct the g^ G SL 2 
that will act on each of the first three qubits. 

Case 1: The set {aj,6j,Cj} spans Hi. Assume {aj, q} is linearly 
independent and 6j = Ajaj+yUjCj. With Aj/Xj ^ we can choose ^ G SL 2 
such that gi(\i<ii) = Vit§ and gi(^iCi) = Vit§ where v\ = det [Aja;|/XjCj]. 
In this situation giipi) = i^(eo + e\). If Ai = then we can choose 
gi G SL 2 such that g;( a «) = u i e o an d 9i{fiiCi) = A similar argument 
holds when ^ = 0. 

Case 2: The set {aj,6j,Cj} does not span Hi. We can assume that 
a i = h = q and choose g± G SL 2 such that gi(cii) = giipi) = gi(ci) = e . 

Now when i = 4 the only difference is that if say {a^, 64} is linearly 
dependent then we cannot simply assume that a 4 = 64, but must take 
into account a scalar factor. 

Figure [3] contains the different ways that the sets {ai, bi, q} can con- 
tain the same vector multiple. It is not surprising that it is more 
complicated than Figure [2] and it is indeed slightly more difficult to 
show that it captures all possibilities. Nevertheless we have: 

Proposition 5.2. Any rank 3 tensor ip G H is SL* oc -equivalent to a 
tensor having one of the patterns (a—g) in Table [3J This pattern is 
uniquely determined by if). 
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Proof. Let if> = a + b + c where 

a = ai (g> a 2 <8> a 3 <g) a 4 , 
b = bi ® b 2 ® b 3 ® 6 4 , 
C = Ci (g) c 2 <g> c 3 <g> c 4 . 

If ■?/> is factorizable, then since it has rank 3, by the classification in [3] 
it must be SL* oc -equivalent to a tensor of the form (g). Assume if) is 
not factorizable. If for each i there are no linear dependencies between 
any two of the factors cij, foj, Cj then if> is clearly in the form (a). If 
there is one linear dependency, by permuting qubits, we may assume 
if> is of the form (b). Assume there are two linear dependencies. If they 
are both between the factors of the same two summands, say a and b, 
then we may assume that b = ai (g a 2 <g> b 3 <g> 6 4 . Now we can use the 
SLi oc -operations as described above to get 

v~ l i\)' = e <S> e <8> e <8> e 

+ e <8> e ® (e + ei) (g> (e + ei) 

+ ei ® ei (g) ei ® e x 

= e (g) e (g) e <E> (2e + ei) 

+ e ® e (g) ei (g) (e + ei) 

+ ei (g ei ® ei ® ei 

which has diagram 



Now we can repeat this process one more time to get a tensor in the 
form 3(e). If there is a dependency between a factor of a and b, and 
another dependency between a factor of b and c, then ^ is in the form 
(c). Assume there are 3 dependencies. If they are between the same 
two summands then if) is no longer a rank 3 tensor. If two are between 
a and b, and one is between b and c, then as before we can bring this 
into the form (e). If we have a dependency between factors of a and 
b, b and c, and between a and c then if) is of the form (d). Assume 
that there are 4 dependencies. Clearly if at least 3 are between two 
summands then if) is not rank 3. Otherwise it is straightforward to see 
that the only possibility is a tensor of the form (/). 

The uniqueness assertion follows by inspection of Table [5j □ 

With this in mind, the normal forms follow. 

Remark 5.3. In one of the cases, the proof below depends on the fol- 
lowing important fact, a special case of [30l Theorem 38.6.1]. Let us 
identify TL with the subspace p of g = sos- Fix if) e p and let O C g be its 
SOg-orbit under the adjoint action. Then each irreducible component 
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of O fl p is a single SO4 x SC^-orbit. Moreover all these components 
have the same dimension. 

Proposition 5.4. A rank 3 tensor if> G TC of the given pattern (see 
Figure^) can be reduced using SLj* oc - operations to the form: 

(a) a|0000) + p{\0) + |1)) 04 + 7|HH), apy ^ 0, 

(6) a(|0000) + |1111)) + |0) <g> (|0) + |1))® 3 ; a + 0, 

(c) a(|0000) + |01) <8> (|0) + (I)) 552 + |1111)), a ^ 0, 

(d) |0000) + |011) ® (|0) + |1)) + |1101), 

(e) a(|0000) + |0011) + |1111)), a ^0, 
(/) |0000) + |0011) + |1110), 

(g) |0000) + |0110) + |1100). 

Proof. Let us outline the procedure in cases (b) and (d). The rest of 
the cases follow from much the same reasoning, although case (b) is 
uniquely non-trivial. Using Figure [3] we see that in the case (b) we may 
assume that 

i/j = a|0000) + /3|0) ® (|0) + |1)) 03 + 7|1111) 

where 0^7 7^ 0. We apply the SL 2 -tranformation aeo — > z/eo and 
7ei — > z/ei to the first qubit to get 

a'(|0000) + |1111)) + (3'\0) ® (|0) + |l)f 3 

where a' = u and /3' = z//?a _1 . 
Let 

^ = a(|0000) + |1111)) + f3\0) ® (|0) + |1))® 3 

with af3 7^ (we rename a! to a and (3 1 to /? for convenience). Choose 
7 G C such that 7 2 + 7 = a(a + /3). If /3 = —a we assume that 7 = —1. 
We claim that if) is SL loc -equivalent to 

= 7 (|0000) + |1111)) + |0) ® (|0) + |1))® 3 . 

If /? = —a then ip = —acf). Since is nilpotent, our claim holds. 

Now assume that a + f3 7^ 0. Choose a continuous function / : 
[0, 1] -> C\{0} such that /(0) = a, /(l) = 7 and /(t) 2 7^ a(a + (3) for 
all £ G [0,1]. Consider the one-parameter tensor family 

xit) = /(t)(|0000> + 11111)) + ^ + P |0) ® (|0) + |i)f 3 . 

It is easy to verify that all matrices R x (t) have the same Jordan struc- 
ture: 

0,J 2 (±Va 2 + «/3),J 3 (0). 
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Hence they all belong to a single GL 8 -orbit (a similarity class) O C M 8 . 
Since q C M 8 is the space of skew-symmetric matrices, O fl g is a single 
Og-orbit, and so it is the union of at most two SOg-orbits. By Remark 
15.31 each irreducible component of O fl p is a single S0 4 x S0 4 -orbit 
and all these components have the same dimension. Since {R x (t)} is 
contained in a single irreducible component of O Hp, and x(0) = ip and 
= 0, we conclude that R^ and i?^ belong to the same SO4 x SO4- 
orbit, and so ip an d belong to the same SLi oc -orbit. This concludes 
the proof of our claim. 

A tensor ip of the form (d) can be reduced to 

a|0000) + /3|011) <8> (|0) + |1)) + 7 |1101). 

By mapping ae^ —>■ ve§ and 7 ei — > ve\ in the third qubit we attain 

a'|0000) + /?'|011) ® (|0) + |1)) + a'|1101) 

and we can check that a' = v and f3' = v(3a~ l . Now by applying the 
SL2-tranformation which sends a'eo —> z/'eo and f3'e\ — > v'e\ in the 
third qubit we obtain 

i/(|0000) + |011) ® (|0) + |1)) + |1101)). 

Since ip is nilpotent, the SLJ" oc and GLj* oc -orbits coincide and we can 
replace v' with 1. □ 

We shall say that the tensors listed in Proposition 15.41 are of type 
3 (a— g), respectively. Note that the invariants L and M vanish on each 
of these tensors. 

For the computation of tensor ranks it is important to know the 
Jordan structure of the matrices for all types of tensors of ranks 
< 3. We shall investigate the tensor ip of type 3 (a) in detail. The 
other cases are easy to analyze and we omit their discussion. By using 
f l3.ll) we find that the characteristic polynomial of R^ is 

t 2 (t 6 + 2(aP + a 7 + /5 7 )t 4 + (a(3 + a 7 + f3j)H 2 + 4(a/3 7 ) 2 ). 

If we let s = t 2 then it is sg(s) where 

g (s) = s 3 + 2(a(3 + cry + /3 7 )s 2 + (a(3 + a 7 + /3 7 ) 2 s + 4(a/5 7 ) 2 . 

The discriminant of g(s) (see (13.21) ) is 

16(«/5 7 ) 2 ((«/? + « 7 + /3 7 ) 3 - 27(«/? 7 ) 2 ). 

If the discriminant does not vanish then R^ is semisimple and belongs 
in family 1 from Table [TJ Now if it vanishes then we must have 

(5.1) (a/3 + a 7 + /5 7 ) 3 -27(a/? 7 ) 2 = 
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and the roots of g(s) are 

(5.2) A 2 = --(a/3 + cry + /i 2 = --(a/3 + a 7 + /3 7 ), 

where the latter is a double root. 

Now we can determine the Jordan structure of R^. Since the eigen- 
value of R^p has multiplicity two, and the matrix rank of R^ is 6 
(there are two rows of zeros) we conclude that there are two lxl 
Jordan blocks of 0. Consider the matrix R^. By permuting rows and 
columns, we see that it is similar to a matrix of the form Pi © © P2 
where the second summand is a 2 x 2 zero matrix and the other two 
summands are 3x3 matrices. It is straightforward to see that Pi and 
P 2 are similar. Now assume that R^ is semisimple, then so is P 2 . In 
particular the matrices Pi — A 2 / 3 and Pi — yU 2 / 3 have ranks 1 and 2 
respectively. But then by evaluating the 2x2 minors we can conclude 
that a — [3 — 7. Conversely one can check that a = (3 = 7 implies 
that R^ is semisimple. We find that P^ then has the Jordan structure: 
0,0, ±A, J 2 (±/i) and it is in family 2, unless a = (3 = 7 in which case 
R^ is semisimple and is again in family 1. 

To summarize, we find that if a tensor of rank 3 is semisimple, then 
it must be of type 3(a). Furthermore, a rank 3 tensor of type (a) is 
semisimple unless ( 15. ip holds and a, /3, 7 are not all equal. 

Table [5] presents the Jordan structures for the different type of ten- 
sors of rank < 3. Note that one has to permute the four qubits in order 
to obtain all possible Jordan structures for a given type. 

TABLE 5. Jordan structure of P^ for tensors of rank < 3 

1(a) 0,0,0,0,J 2 (0),J 2 (0). 

2(a) 0,0,0,0,±zv^,±iv^. 
2(6) 0,0,J 3 (0),J 3 (0). 

2(c) 0,0,0,0,0,J 3 (0) and J 2 (0), J 2 (0), J 2 (0), J 2 (0). 

3(a) Discussed above: Semisimple or 0,0, ±zA, J 2 (±i/x). 

3(6) If a = -1 then 0, J 7 (0) else 0, J 2 {±Wa 2 + a), J 3 (0). 

3(c) ±ia,±ia,Q, J 3 (0) and J 2 (0), J 2 (0), J 2 {±ia). 

3(d) J 3 (0),J 5 (0). 

3(e) ±ia, ±ia, J 2 (0), J 2 (0) and 0,0,0,0, J 2 (±ia). 

3(f) 0,0,0, J 5 (0) and J 4 (0), J 4 (0). 

3(^) 0,J 2 (0),J 2 (0),J 3 (0). 
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6. Determining the tensor ranks 

Here we will compute the ranks of tensors if) in each of the nine 
families listed in Theorem 13.61 For any subsequence {a^} of 1,2,3,4 
we set H ax ... ak = <S>^ =1 Ti. ai . For if> = e <E> to + ei (g> t\, where {t , *i} is 
linearly independent, we form the linear transformation : C 2 — > H234 
sending (x, y) — > The image of C 2 \{0} under is a projective 

line in the projective space P(7i234)- For a pure state G H234 we 
will denote by det0 the Cayley hyperdeterminant as described in [3]. 
Explicitly we have 

det = (tr A tr B - tr AB) 2 - 4 det A det 5 

where 



A = 


0000 


0001 


; B = 


0100 


0101 




_0O1O 


001 1_ 




0110 


0111 



Note that 4| det0| is the residual entanglement (also known as the 3- 
tangle) of the pure state described in [6]. 

Lemma 6.1. If (xi,yi) and (x2,y 2 ) are linearly independent then we 
have that rank-0 < rankT^(;ri, yi) + rankT^(x2, y 2 )- 

Proof. Since {(x 1 , yi), (x%, 2/2)} is linearly independent, so is {xie + 
^2ei,2/ie +y 2 ei}. The tensor ■0' = (x 1 eo+x 2 ei)®t o + (y 1 e o + y 2 e 1 )0t 1 is 
in the same GLi oc -orbit as -0 so rank?/*' = rank-0. But if>' = eo® (xit + 
yih) + ei <g> (x 2 to + Z/2*i) so rank-0' < rankT^,(xi, yi) + rank T,/, (x 2 , 7/2) 
and the result follows. □ 

Table lists, in order, the invariants H,L,M,N, D and T for each 
of the 17 families from Table [TJ It will be helpful to refer to the list 
of invariants of the families as we proceed. We will often use the fact 
that if either L or M does not vanish on ip then if) ^ £3 and so if} must 
have rank 4 (see Proposition 13.11) . It is easy to determine whether a 
state if) is factorizable. For instance, we have if> — <fi ® % with e Jtvi 
and x ^34 iff rank?/> = 1. Similarly we have if) = <g> x with e Hi 
and x £ ^234 iff the 2x8 matrix [V'ijfcz] has rank 1. In general, one 
has first to permute the qubits. 

We shall now consider separately each of the nine families mentioned 
in Theorem 13.61 

Family 1. We may permute the diagonal entries a,b,c,d of R (see Table 
[I]) and replace them by ±a,±b,±c,±d without changing the SLj* oc -orbit 
of if). If L 7^ or M 7^ then rank-0 = 4. From now on we may assume 
that a + b + c = d = (see the expressions for L and M in Table [6]) . 
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Table 6. The SLi oc -invariants H, L, M, N, D, F 

1. \{a 2 + b 2 + c 2 + d 2 ) ; abed ; ±(4(ad - be) 2 - (a 2 -b 2 -c 2 + d 2 ) 2 ); 
-±(4(ac + bd) 2 -(a 2 -b 2 + c 2 -d 2 ) 2 ) ; -\(ad-bc)(ab-cd)(ac-bd); 
-^(2(a 6 + b e + c 6 + d 6 ) - (a 2 + b 2 + c 2 + d 2 )(a 4 + o 4 + c 4 + d 4 ) + 
18(aW + b 2 c 2 d 2 + c 2 d 2 a 2 + d 2 a 2 b 2 )) 

2. \{a 2 + b 2 + 2c 2 ) ; abc 2 )~(a- 6) 2 ((a + bf - 4c 2 ); i(a + o) 2 ((a - 
o) 2 - 4c 2 ) ; 

ic 2 (a - o) 2 (c 2 - 06) ; ^ ( (a 2 + o 2 ) ( (a 2 - o 2 ) 2 + 16c 4 ) - 2c 2 (a 4 - 18a 2 6 2 + 
& 4 )) 

3. i(a 2 + b 2 ) ; ; -^(a 2 - o 2 ) 2 ; ^(a 2 - b 2 ) 2 ; ; ^(a 2 - b 2 )(a 4 - o 4 ) 

4. i(a 2 + 6 2 ) ; ; -^(a 2 - 6 2 ) 2 ; ^(a 2 - o 2 ) 2 ; ; ^(a 2 - o 2 )(a 4 - o 4 ) 

5. a 2 + b 2 ■ a 2 b 2 ■ ; -a 2 b 2 ; ; a 2 b 2 {a 2 + b 2 ) 

6. ±(a 2 + 36 2 ) ; ab 3 ; ^(6-a) 3 (a + 36); ±(a -3b)(a + b) 3 ; ±o 3 (o-a) 3 ; 
^((a 2 -o 2 ) 3 + 16o 4 (3a 2 + o 2 )) 

7 I a 2 . Q . _J_ a 4 . J_ fl 4 . Q . J_ 6 

2 ' ' 16 '16 ' ' 32 

8. ~^2^ 7 7 -^g a , a , j 32~^ 

9. 2a 2 ; a 4 ; ; -a 4 ; ; 2a 6 

10. a 2 ; ; ; ; ; 

11. a 2 ; ; ; ; ; 

* Families 12-17 are nilpotent so all invariants are 0. 



Then we have 

Aip = a(e + ex) ® (e + e x ) ® (e Q - e x ) ® (e - ei) 
+ a(e - ei) ® (e - ei) ® (e + e x ) ® (e + ei) 
— 4c(e ® ei (g> ei <8> e + c\ <S> e <8> e <S> ei). 

If a&c = 0, say c = 0, then rank?/; = 2. Note that after normalization, 
this ip represents the generalized GHZ state. Now we may assume that 
abc 7^ 0. Let <fi be the tensor of type 3(a) (see Proposition 15. 4p . We 
shall choose the scalars a, /3, 7 to satisfy the two equations 

(6.1) 2(a/5 + /?7 + 7a) = a 2 + 6 2 + c 2 , 4(a/?7) 2 = a 2 6 2 c 2 . 

and to ensure that is semisimple. If a, b, c are not distinct, say a = b, 
we can take a = (3 = 7 = a. For semisimplicity of <fr m this case see 
the end of the previous section. 

From now on we assume that a, 6, and c are distinct. If a 2 +6 2 +c 2 = 
we take (3 = a( and 7 = a( 2 ; where ( = e 2 ™/ 3 , and we choose a such 
that —2a 3 = abc. Otherwise we take (3 = —a, 7 = and choose a 
such that —2a 2 = a 2 + b 2 + c 2 . The equations ( 16. II) imply that and 
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Rtj, have the same characteristic polynomial. Since the nonzero eigen- 
values of R^p, i.e. ±ia,±ib,±ic, are distinct, one can verify easily that 
R^ is semisimple. As L and M vanish on ip and </> and the equations 
( 16.11) show that H and D also agree on ip and 0, Theorem 13.31 implies 
that ij} and <fi are in the same SLi oc -orbit. Hence rank-0 = rank(/> = 3. 

Family 2. If ip ^ S3 then rank-?/; = 4. Otherwise we must have 
L = M = 0, i.e., afec = (a - 6) ((a + b) 2 - Ac 2 ) = 0. As we may 
switch a and b and multiply a, b, c by ±1, there are only three cases to 
consider: 

(i) a = b, c = 0; 
(iij a = b = 0, c ^ 0; 
(iii) 6 = 0, a = -2c ^ 0. 

In case (i) we have 

Aip = (e - e x ) (g (e + ei) <g> (e + ei) <g> ((a - 2z)e + (a + 2z)ei) 
+ a(e - ei) ® (e - e x ) <g (e - ei) <g (e - ei) 
+ 2ae x <g> (e + ei) (g (e + ei) <g> (e + ei), 

and in case (ii) 

Ai[) = c(e + ei) <g> (e - ei) g) (e + e x ) (g> (e - e x ) 
+ c(e - ei) g) (e + d) g) (e - ei) g) (e + ei) 
- 2z(e - ei) g) (e + ei) g) (e + ei) ® (e - ei). 

Clearly if a = in case (i) then rank?/' = 1. Otherwise in cases (i) 
and (ii) it is easy to verify that ip ^ S2 and so rank?/' = 3. Now we 
consider case (iii). Let <fi G H be of type 3(a) with a — |, /3 = zca/3 
and 7 = — /?. Then the matrices and Pi?,/,, where P is the diagonal 
matrix with (—1, 1, 1, 1) diagonal entries, are symmetric and have the 
same Jordan structure: 0, —2c, Ja(c). Hence they are orthogonally 
similar. This shows that R^ and R^ are in the same O4 x 04-orbit. 
Hence ip and <p are in the same SL* oc -orbit and so rank?/) = rank</> = 3. 

Family 3. If a 2 7^ b 2 then ip ^ 5*3 and rank?/> = 4. Since we can 
interchange a and 6 and replace them by ±0 and ±6, we may assume 
that a = b. If a = then 

2-0 = (e + ei) (g (e - e x ) (g (ei (g e - e g) e x ) 
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and rank-0 = 2. If a 7^ then if) ^ 62 and 

4aip = (e + ei)® ((a 2 + l)e + (a 2 -l)ei)® (e + ei)®(e + ei) 

- (eo + ei) <g> (e - ei) ® ((1 + a)e + (1 - a)ei) 

®((1 - a)e + (1 + a)ei) 

- 2a 2 d <g) (e - ei) <g> (e - ei) <g> (e - ei), 

and so rank?/' = 3. 



Family 6. If if) (jL S3 then rank-?/; = 4. Otherwise L = M = and 
by using Table El we have a = b = 0. Then if) is nilpotent and it is easy 
to verify that if> G §2- Since ^ is not factorizable, it cannot have rank 
1 or 2 (see Proposition [5J]). The matrix R^ has Jordan structure 0, 0, 
^3(0), ^3(0) (see Tabled]). Since this is absent from the rank 3 section 
of Table 0, we infer that rank-?/; 7^ 3. Thus rank?/> = 4. 

Family 9. If a 7^ then if) ^ S3 and so rank?/) = 4. If a = then 

if) = 2z(|1110) - |0010) + |1001>) 
and if) $l §2 so rank-?/' = 3. 

Family 10. Permute qubits 1 and 2. (The effect on the matrix is just 
to change its (3, 3)-entry from 1 to —1.) If a 7^ then rankT(l, 1) = 1 
and rankT(l, 0) = 2. Since if) ^ 5*2 we have that rank?/> = 3. If a = 
then if> — (ei — e ) <g> <fi. Since det = and <f> is not factorizable, we 
have rank?/> = rankc/> = 3. 

Family 12. We check that if) ^ S~2 since the 1,1 minor of if) is nonzero. 
We have that rankT(l, 0) = 1 and rankT(l, 1) = 2. Hence rank-0 = 3. 
A computation gives 

(\-i)\plif) = (1 — i)\/2(eo - ei) (g) ei ® (e — iei) (g) (— ie + ei) 

- ei <g> (eo - i\/2ei) <g> (e + /3ei) <g> (e Q + aei) 

- ei ® (e + i\/2e x ) ® (-e + ae x ) <g> (e - /3ei), 

where a = \/2 + 1 and (3 = a/2 — 1 . 

Family 14- It is again straightforward to verify that ip ^ §2- Then using 
that rankT(i, 1) = 1 and rankT(l, 0) = 2, we obtain that rank-?/> = 3. 



Family 16. In this case if) = (eo + ei) (g> 0. Since det</> ^ we have 
rank-0 = rank</> = 2. 
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After all these computations it is worthwhile observing that S3 con- 
tains only one SLj* oc -orbit of rank 4 tensors. This exceptional orbit is 
the unique nilpotent orbit of family 6 (with a = b = 0). It is the orbit 
of the generalized W state given by (14. 2 p and it is contained in S 2 - In 
particular we have S2 7^ S2 and S3 7^ S3. 

7. Tensor rank algorithm 

By using the results of the previous section, we can now construct 
a simple algorithm for computing the tensor rank. We have explained 
in the previous section how to test a state if) for factorization. If one 
of the factors is from a single Tik, we may use density matrices. For 
a state if) G 7i let p = \ip)(if>\ be its density matrix. Denote by pk its 
reduced density matrix obtained by tracing out all qubits but the k-th 
one (for the definition of the density matrices and partial trace see e.g. 
[27]). Then if) factorizes, with one of the factors in Tik, iff the matrix 
Pk has rank 1, so we let be the matrix rank of pk- With an abuse 
of notation, we let rk be the rank of the corresponding pk for 3-qubit 
tensors as well. 

We now give our algorithm for computing the tensor rank of an 
arbitrary state if> G H- The algorithm uses another procedure which 
computes the tensor rank of 3-qubit states, which can be deduced from 
[3]. It should be understood that the algorithms halt as soon as the 
rank is returned. Recall the definition of the hyperdeterminant det if) 
for 3-qubit pure states if) given in the previous section. 

3— Qubit Tensor Rank Algorithm. 

Input : A nonzero tensor if> G 7^123 
Output : The tensor rank of if) 

1 If det if) is nonzero then return 2. 

2 Compute the ranks of pk for k G {1,2,3}. 

3 If Tk = 1 for at least two different k then return 1. 

4 If some = 1 then return 2 . 

5 Return 3 . 

With this we may compute ranks of 4-qubit tensors. 

4- Qubit Tensor Rank Algorithm. 

Input : A nonzero tensor if) G H 
Output : The tensor rank of if) 

1 If L(if)) or M(if>) is nonzero then return 4. 
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2 If at least one of the forty-eight 3x3 minors of the 
matrices ip , ip' , ip" is nonzero then return 3. 

3 Compute the ranks r& of . 

4 If say n = 1, then ip — V\ ® (p with V\ G T~L\ and <p G 7^234 ) 
and return rank0. 

5 Now all r^ = 2. If ip is nilpotent, i.e., i?^ is nilpotent, 
then return 4. 

6 Return 2 . 

Let us show that the algorithm is correct. It may be helpful to look 
at Figure H] where some of the sets we use below are exhibited. Step 1 
is clear. In order to justify step 2, it suffices to verify that 

ip G S 3 \S 2 =>- rank-?/* = 3. 

This follows from the case-by-case analysis of the previous section. 

After reaching Step 3, we have ip G S2. Consequently, the families 
7,8,9,12,13,14 and 15 are ruled out, i.e., ip does not belong to any of 
them. Indeed it is easy to verify that none of these families meets §2- 

Steps 3 and 4 are also clear. 

After reaching step 5 our ip is non-factorizable and we can rule out 
the families 2,3,4,5,10,11,16 and 17. Indeed the families 16 and 17 are 
factorizable and the orbits in the families 2,3,4,5,10 and 11 which are 
contained in S2 are also factorizable. Hence ip belongs to the family 
1 or 6. If it is nilpotent, it is in family 6 and has rank 4. Otherwise 
it is in family 1 and the detailed analysis of this case in the previous 
section shows that the rank of ip is 2. 

Figure 0] describes the structure of Ti. with respect to tensor ranks. 
Each vertex represents a Zariski closed set and it is ordered by inclusion 
as one progresses to the top vertex. All sets on or below the horizontal 
line satisfy the equation H = and consist of nilpotent orbits. However 
there exist nilpotent orbits not contained in S 2 - The numbers between 
pairs of adjacent vertices indicate the rank of the tensors that are in 
the set theoretic difference between the higher and lower vertices. The 
numbers on the far left indicate the dimension of the corresponding 
affine variety. 

8. Conclusion 

In this paper we have investigated the SLOCC classification of pure 
states of four qubits first described by Verstraete et al. in [31]. The 
families of representatives provided in that paper were accurate except 
for some possible misprints in the family L 6 3 . However their claim of 
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uniqueness in Theorem 1 is not true and the subsequent proof was not 
easy to follow. 

We have provided a more general version of that theorem in our 
Theorem 12.81 We presented this theorem within the framework of 
orthogonal representations of a certain quiver Q. We felt that this 
approach lead to a simpler proof of Theorem 12.81 We also observed 
that this theorem can be deduced (with some additional work) from the 
theory of symmetric quivers as presented in a recent paper of Derksen 
and Weyman [7]. 

We found it beneficial to embed the 4-qubit Hilbert space 7i into 
the Lie algebra g of the complex orthogonal group 0s- This naturally 
lead to the notion of semisimple and nilpotent states. The semisimple 
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states are dense in H, while the nilpotent ones comprise only finitely 
many SLi oc -orbits. The subgroup 4 x 4 C 8 acts naturally on H. 

We have also provided a more complete description of the behaviour 
of the 4 x 04-orbits under permutations of qubits. This was ad- 
dressed in [31] as well, but an important characteristic was not stressed. 
Namely that the action of a permutation of qubits on a family of orbits 
may map some orbits into a different family while at the same time 
map others back into itself. So in general, a permutation of qubits 
does not induce a permutation of families of orbits as naivete would 
have one think. 

The problem of showing that two states 0, ip G H are not SLj* oc - 
equivalent appears in a number of recent papers [2J [33]. We derived 
polynomial invariants for SLj* oc and show that two semisimple states 
are SL* oc -equivalent iff they agree on the invariants. In 1852, these 
same invariants were considered by Schlafli who noted their invariance 
under permutations of indices. The general case is somewhat more 
complicated as it requires computing the Jordan structure of associated 
matrices R, and possibly the use of a6-diagrams. 

The other focus of this paper was to ultimately develop an algorithm 
that would calculate the tensor rank for a nonzero tensor in 7i. This 
was accomplished by a thorough examination of each of the families in 
the SLj* oc classification. To carry out this analysis we used the results 
of Brylinski in [3]. In particular it was essential to know that the 
maximum rank of a tensor in 7i is 4 and that the polynomial SLi oc 
invariants L and M define the Zariski closure of the tensors of rank 
< 3. We found that another set of 48 equations define the Zariski 
closure of the tensors of rank < 2 based on the speculation by Brylinski 
in [3]. It was also fortunate that the tensors of rank < 3 admitted a 
simple classification which allowed us to deduce the ranks of some other 
tensors in certain cases. We then were able to construct the algorithm 
which is pleasantly simple compared to the machinery involved in the 
analysis mentioned above. 

The authors of [31] claim to have solved the problem of equivalence 
of two states under the group Ui oc of local unitary operations. For that 
purpose they propose a normal form via a two step procedure. Appar- 
ently they failed to observe that the second step of their procedure may 
undo the beneficial effect of the first step. It is unclear how their pro- 
posed normal form is actually defined. Since Ui oc is a compact group, 
this equivalence problem can be solved by computing the algebra of real 
polynomial invariants 7i — > K for Ui oc . Indeed it is known in general 
that these invariants separate Ui oc -orbits. However, a set of generators 
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for this algebra of invariants has not been computed so far although 
its Poincare series has been computed independently in [32] and [2"U] . 
The problem of local unitary equivalence is also considered in [22] but 
the results are far from conclusive. Hence this problem remains open. 

Appendix A 

The following table gives the representatives ip °f the nine fami- 
lies of SL* oc -orbits of non-normalized pure states of four qubits. They 
are derived from the corresponding i?-matrices given in Tabled] by the 
transformation — > ip — T~ 1 R^T. We point out that our .R-matrices 
for these families are chosen to be as simple as possible, consequently 
the expressions for the corresponding states ip are not. The represen- 
tations given in [3_T] are in some cases shorter than our representations, 
e.g. for family 16 their representative is 

|0)®(|000) + |111)). 

We recall that SL* oc -orbits originating from different families are 
necessarily distinct, and that two states in the same family may be in 
the same SLf -orbit. 
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Table 7. Representatives ip of SL* oc -orbits 



1. ^(|0000) + |1111)) + ^(|0011) + |1100)) 
+*±£(|0101) + |1010)) + ^flOHO) + |1001>) 

2. ^^(|0000) + |1111>) + 2^±*(J0011) + |1100>) 
+%ti(|0101) + |1010)) + ^p(|0110) + |1001)) 
+|(|oooi) + join) + 1 1000) + 1 1110) 

-|0010) - |0100) - |1011) - |1101>) 

3. § (|0000) + 1 11 11> + |0011) + 1100)) + *±!(|0101) + |1010>) 
+*§±(|0110) + |1001)) + |(|1101) + |0010) - |0001) - |1110)) 

6. ^(|0000) + 1 11 11>) + &Q0101) + |1010)) + i(|1001) - |0110>) 
+ ^(|0011) + |1100>) + Klooio) + |0100) + 1 1011) + 1 1101) 
— 10001> - |0111) - 1 1000) - |1110>) 

9. a(|0000) + |0101) + |1010) + 1 11 11>) 
— 2f(|0100> - |1001) - 1 1110>) 

10. s±i([0000> + 1 1 11 1> + |0011) + 1100)) + ^=|±1(|0101) + |1010>) 
+^=|=1(|0110) + |1001>) + i±i(|1101> + |0010)) 
+i§±(|0001) + |1110>) - |(|0100) + |0111) + |1000) + |1011>) 

12. (|0101) - |0110) + |1100) + 1 11 11>) + 1)(|1001) + |1010>) 
— z(|0100> + |0111) + |1101) - |1110>) 

14. ^±i(|0000) + |1111) - |0010) - |1101)) 
+i§l(|0001) + 1 1110) - |0011) - 1100)) 
+ |(|0100) + |1001) + |1010) + |0111)) 
+±=2(11000) + |0101) + |0110) + |1011>) 

16. |(|0) + |1)) <g> (|000) + |011) + |100) + I 111) 
+i(|001) + |010) - 1 101> - 1 110))) 



Appendix B 

To simplify the notation, we shall use the same symbols to de- 
note the polynomials in A or A* and their restrictions to the Car- 
tan subspace a. As a byproduct of our construction of generators of 
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the algebra A*, we have obtained a nice set of generators of the al- 
gebra of polynomial invariants of the Weyl group of type _F 4 . Their 
degrees are, of course, 2, 6, 8 and 12. The invariant of degree 12 
has the factorization II = (L - M)(M — N){N - L) and the one 
of degree 8 is a sum of three squares E = L 2 + M 2 + iV 2 where 
L = abed, 

M = ±(A(ad - be) 2 - (a 2 - b 2 - c 2 + d 2 ) 2 ) 

= -±((a + b) 2 -(c + d) 2 )((a - bf - (c - d) 2 ), 
N = -±(4{ac + bd) 2 -{a 2 -b 2 + c 2 -d 2 ) 2 ) 

= i((a + b) 2 - (c - d) 2 ){{a - b) 2 - (c + d) 2 ). 
The known sets of generators which we could find in the literature 
[T5l GU [28] do not share these special features. 

Table 8. Generators H, T, E, II of invariants of the 
Weyl group of type F4 



2H = a 2 + b 2 + c 2 + d 2 

2 5 r = 2(a 6 + 6 6 + c 6 + rf 6 ) - (a 2 + b 2 + c 2 + d 2 )(a* + b 4 + c 4 + d 4 ) 

+18(aW + b 2 c 2 d 2 + c 2 d 2 a 2 + d 2 a 2 b 2 ) ' 
2 7 E = 2 8 a 2 b 2 c 2 d 2 + {A{ad - be) 2 - (a 2 - b 2 - c 2 + d 2 ) 2 ) 

(4(ac + bd) 2 - (a 2 - b 2 + c 2 - d 2 ) 2 ) 
2 12 n = (Wabcd - A{ad - be) 2 + (a 2 - b 2 - c 2 + d 2 ) 2 ) 

{A(ad - be) 2 - (a 2 - b 2 - c 2 + d 2 ) 2 

+4(ac + bd) 2 - (a 2 - b 2 + c 2 - d 2 ) 2 ) 

(-4(ac + M) 2 + (a 2 - 6 2 + c 2 - d 2 ) 2 - lQabcd) 

The generators 7 2 , hi h an d /12 found in relate to our generators 
as follows: 
h = 12H, 
h = 72H 3 - 96r, 
Is = 2UH 4 - 832r# + 320E, 

7 12 = 4104# 6 - 24096F 3 r + 17440F 2 E + 3904r 2 - 3840n. 
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